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Abstract A review of some recent results on the dynamical theory of the Yang- 
Baxter maps (also known as set-theoretical solutions to the quantum Yang-Baxter 
equation) is given. The central question is the integrability of the transfer dynamics. 
The relations with matrix factorisations, matrix KdV solitons, Poisson Lie groups, ge- 
ometric crystals and tropical combinatorics are discussed and demonstrated on several 
concrete examples. 



Introduction 

The story of Yang-Baxter maps is usually traced back to Drinfeld [1] , who sug- 
gested to study the set-theoretical solutions to the quantum Yang-Baxter equa- 
tion, but the first interesting examples were constructed already by Sklyanin [2] . 
Drinfeld's question was immediately appreciated by Weinstein and Xu [3] who 
developed an approach to this problem based on the theory of Poisson Lie groups 
and symplectic groupoids. Etingof, Schedler and Soloviev [4,5], partly inspired 
by Hietarinta's work [6], took more algebraic point of view (see also [7-12]). 

Very interesting examples of such solutions have appeared in the theory of 
geometric crystals (Berenstein-Kazhdan [13] and Etingof [14]). This opened re- 
markable links with the theory of soliton cellular automata (" box-ball" systems) 
and tropical combinatorics extensively developed by the Japanese school [15-23]. 

The study of the dynamical aspects of this problem was initiated in the pa- 
per [24], where the term "Yang-Baxter map" was introduced. My motivation 
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came from the soliton theory, in particular from the theory of matrix KdV equa- 
tion [25]. I was puzzled by what kind of information one can learn from the in- 
teraction of two matrix solitons and what is the role of the Yang-Baxter relation 
here. Analogy with the theory of solvable models in statistical mechanics and 
quantum inverse scattering method [28-32], where the quantum Yang-Baxter 
equation has its origin, naturally led to the notion of the transfer dynamics and 
to the relations with the theory of integrable systems. The main idea of [24] 
was that the investigation of the transfer dynamics should reveal the structures 
encapsulated in a given Yang-Baxter map. A very good example is the Adler 
map (see below) from the theory of the dressing chain [33,34]. The correspond- 
ing transfer dynamics is nothing else but the discrete version of the finite-gap 
KdV dynamics (Novikov equations [35,36]) with its rich algebro-geometric and 
symplectic structures. 

In this extended version of the lectures given at Kyoto meeting " Combina- 
torial Aspects of Integrable Systems" in July 2004 I am going to review the 
dynamical theory of the Yang-Baxter maps and illustrate it on several concrete 
examples following mainly to [24-27]. 

We start with the general discussion of the Yang-Baxter maps and related 
transfer dynamics. The main property of the transfer maps, which they share 
with the transfer matrices, is the pairwise commutativity. If the maps are 
polynomial or rational the commutativity relation is quite strong and sometimes 
implies the solvability of the corresponding dynamical system (see [37] for the 
discussion of some known results in this direction). Although it is not clear 
what is the transfer dynamics in general, our main motive is that it should be 
"integrable" and one can think of it also as a set-theoretical analogue of the 
finite-gap theory. 

We illustrate this by two constructions of the Yang-Baxter equation com- 
ing from the theory of integrable systems. The first one is related to matrix 
factorisations and QR-type procedure, which is a discrete analogue of the Lax 
representation. We explain how the Lax matrices can be read off the map itself 
following the procedure suggested in [26] . 

The second class is given by the interaction of solitons with non-trivial inter- 
nal parameters (e.g. matrix KdV solitons). It is a well-known phenomenon in 
soliton theory (see e.g. [38-42]) that the interaction of n solitons is completely 
determined by their pairwise interactions (so there are no multiparticle effects) . 
The fact that the final result is independent of the order of collisions means that 
the map determining the interaction of two solitons satisfies the Yang-Baxter 
relation. 

We briefly consider also the examples of the Yang-Baxter maps related to 
geometric crystals and tropical combinatorics [14, 19, 22]. Then we discuss 
the Hamiltonian theory of the Yang-Baxter maps following recent work by 
Reshetikhin and the author [27] partly based on some general considerations 
from [3]. The main observation is that if the Lax matrices form symplectic 
leaves in some Poisson Lie groups, then the corresponding transfer dynamics 
is symplectic. In the last section we discuss the classification problem of the 
Yang-Baxter maps and its relation with the classification of the quadrirational 
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maps done recently by Adler, Bobenko and Suris [43]. 

Yang-Baxter maps and transfer dynamics 

The original quantum Yang-Baxter equation is the following relation on a linear 
operator R:V ®V -s- F ®V : 

R12R13R23 = R23R13R12, (1) 

where Rij is acting in i-th and j'-th components of the tensor product V ®V ®V 
(see e.g. [32]). 

Following Drinfeld's suggestion [1] consider the following set-theoretical ver- 
sion of this equation. 

Let X be any set, R : X x X — > X x X be a map from its square 

into itself. Let R i3 ■ : X n -> X n , X n = X x X x x X be the maps 

which acts as R on i-th and j-th factors and identically on the others. More 
precisely, if R{x,y) = (f(x,y),g(x,y)),x,y e X then R lj (xi,x 2 , ■ ■ ■ ,x n ) = 
(xi,x 2 , . . . ,Xi-i, f(xi,Xj),x i+ i, . . . ,Xj-i, g(xi,Xj),Xj + i, . . . ,x n ) for i < j and 
(xi,X2, ■ ■ ■ ,Xj-i,g(xi, Xj),x j+ i, . . . ,Xi-i, f(xi,Xj),Xi+i, ...,x n ) otherwise. In 
particular for n = 2 R2i(x,y) — (g{y,x), f(y,x)). If P : X 2 — > X 2 is the per- 
mutation of x and y: P(x,y) = (y,x), then obviously we have 

R21 = PRP. 

Following [24] we will call R Yang-Baxter map if it satisfies the Yang-Baxter 
relation |T]) considered as the equality of the maps of X x X x X into itself. If 
additionally R satisfies the relation 

R 2X R = Id, (2) 

we will call it reversible Yang-Baxter map. 

The condition ^ is usually called the unitarity condition, but in our case 
the term reversibility is more appropriate since in dynamical systems termi- 
nology this condition means that the map R is reversible with respect to the 
permutation P. 

Note that if we consider the linear space V = C x spanned by the set X, then 
any Yang-Baxter map R induces a linear operator in V ® V which satisfies the 
quantum Yang-Baxter equation in the usual sense. Therefore we indeed have 
a very special class of solutions to this equation. If X is a finite set and R is 
a bijection then we have permutation-type solutions discussed in [6]. However 
this point of view on the Yang-Baxter maps seems to be artificial because it 
does not reflect the nature of the maps (like birationality). 

Remark. We should warn the reader that some authors (see e.g. [7, 12, 14]) 
consider instead of the map R a closely related map S = PR, where P is the 
permutation map. The corresponding maps satisfies instead of the Yang-Baxter 
relation the braid relation 

S12S23S12 = S23S12S23 (3) 
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and instead of reversibility the involutivity condition 

S 2 = Id (4) 

One should be careful here since the groups generated by the Rij and S%j maps 
could be different. For example, the group generated by R in general is Z while 
for S it is Zi 2 . This is an important difference from dynamical point of view (see 
below). 

We can represent the relations JIJ [2]) in the standard diagrammatic way as 
shown on the Fig. 1 and 2. 




Figure 1: Yang-Baxter relation 




i ► 



Figure 2: Reversibility 

One should mention that there exist another (dual) way to represent the 
Yang-Baxter relation, which shows the relation with the so-called 3D consistency 
condition in the theory of integrable discrete equations on quad-graphs [44- 
46]. In this representation the fields (elements of X) are assigned to the edges 
of elementary quadrilaterals rather than to the vertices. The corresponding 
"cubic" representation of the Yang-Baxter relation is illustrated in the Fig. [3l 

The left-hand side of ([T]) corresponds to the chain of maps along the three 
rear faces of the cube on Fig. [3J 

#12 : {x,y) (£2,2/1), R13 ■ (x 2 ,z) i-> (x 23 ,zi), #23 : (2/1,21) ^ (2/13,212), 

while its right-hand side corresponds to the chain of the maps along the three 
front faces of the cube: 

#23 : (y,z) )->• (2/3,22), R 13 : (x,z 2 ) ^ (x 3 ,z 12 ), R 12 : (x 3 , y 3 ) >->• (x 2 3, Viz)- 

So, ([l]) assures that two ways of obtaining {x 23 , j/13, z\ 2 ) from the initial data 
(x, y, z) lead to the same results. 
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One can introduce also a more general parameter-dependent Yang-Baxter 
equation as the relation 

#12(Al, A 2 )-Ri3(Ai, A 3 )i?23(^2, A3) = i?23(A 2 , A3)i?i3(Ai, A3)i?i2(Ai, A 2 ) (5) 

and the corresponding unitarity (reversibility) condition as 

R 21 (p,X)R(X^)=Id. (6) 

Although it can be reduced to the usual case by considering X = X x C and 
R(x,\;y 7 [i) — R(X, n)(x, y) often it is useful to keep the parameter separately 
(see the examples below). 

In the theory of the quantum Yang-Baxter equation the following transfer 
matrices 

t (n) . y®n y®n play a f un d a mental role. They are defined as the 
trace of the monodromy matrix 

£(«) — trv„Rt3nRon-l ■ ■ ■ Roi (7) 

with respect to the additional space Vq. If the solution of the quantum Yang- 
Baxter equation depends on an additional spectral parameter A then the corre- 
sponding transfer matrices commute: 

i(">(A)t< n >(/i) = i (n) (M)t (n) (A), (8) 

which is the crucial fact in that theory (see [29,30,32]). 

In the general set-theoretical case we have a problem with the trace operation 
(see although the remark at the end of this section), so one should look for an 
alternative. A possible candidate is given by the operators first considered by 
C.N Yang in the classical paper [28]. Their relation with transfer matrices are 
nicely explained in Gaudin's book [31] (see especially Chapter 10, sections 2 and 
3). A related constructions were used by Frenkel and Reshctikhin in the theory 
of q-KZ equation [47] and by Fomin and Kirillov in the theory of Schubert 
polynomials [48] (see also [49,50]). 




Figure 3: "Cubic" representation of the Yang-Baxter relation 



5 



Following [24] we introduce the transfer maps , n as the maps 

of X n into itself defined by the following formulas: 



T, 



(n) 



ii-\-n—2 • 



R 



(9) 



where the indices are considered modulo n with the agreement that we are using 
n rather than 0. In particular T[ n ^ = i?i„i?i„_i . . . R\2- 

Theorem [24]. For any reversible Yang-Baxter map R the transfer maps 



T^ n \ i = 1, . . . ,n commute with each other: 



rp{n)j,(n) 



(10) 



and satisfy the property 



T[ n) Ti n) . . .rw = Id. (11) 

Conversly, suppose that the maps T. determined by the formula |P|) commute 
and satisfy the relation HI}) for any n > 2 then R is a reversible Yang-Baxter 
map. 

Proof of the first part is very similar to the proof of the commutativity of 
the transfer matrices and follows from the consideration of the corresponding 
diagrams representing the products T> n 'T^ and T^ n 'T> n ' respectively (cf [32]): 




i+ 1 




i + i 



j + i 



Figure 4: Commutativity of the transfer maps 



One can easily check that the second diagram is a result of the several 
operations presented above on figures 1 and 2 applied to the first diagram. The 
second identity also follows from similar consideration (see fig. below where the 
case n = 3 is presented). 

To prove the second part it is actually enough to consider only the cases n = 2 

(2) (2) 

and n = 3. When n — 2 we have two maps T± = R and T 2 = R21, so the 
relation (llip becomes simply the reversibility condition ([2]) . The commutativity 



G 




Figure 5: Product of transfer maps 



condition is trivial in this case. For n = 3 then we have three transfer maps: 
if 5 = R 13 R 12 ,T^ } = i? 2 ii? 23 ,r 3 (3) = iiaaiZsi. The product T^T^T^ = 
R13R12R21R23R32R31 is obviously an identity because R12R21 = R23R32 = 
R 13 R 31 = Id. Commutativity T 1 (3) T 2 (3) = T 2 (3) T 1 (3) means that R13R12R21R23 = 
R21R23R13R12 which is equivalent to Yang-Baxter equation |T]) modulo the 
reversibility relation which was already shown. 

The properties of the transfer maps can also be explained in terms of the 
corresponding extended affine Weyl group A„_ x . As we have already mentioned 
above the maps 

Si = Pii+iRu+i, 

where is the permutation of z-th and j-th factors in X n , satisfy the relations 
SiSi+iSi — Si+iSiSi+i, i = l,...,n — l (12) 

and 

S? = Id (13) 

which are the defining relations of the affine Weyl group A^_ v However (see 
the warning after the definition of the Yang-Baxter maps) the transfer maps 
do not belong to this group but do belong to its extension A^li generated by 
Si,...,S n and the cyclic permutation u = Pi n Pin-i ■ ■ ■ P12, oj n — Id (known 
as extended affine Weyl group, see [51]). Moreover one can easily check that they 
corresponds to the translations in this groups and thus must commute (cf. [50] ) . 

Note that if we consider only the set of Si with i = 1 , . . . , n — 1 (without 
S n ) then they generate the (finite) permutation group S„. In that case as it was 
shown by Etingof et al [4] the corresponding dynamics is equivalent to the usual 
action of S n by permutations (and thus is trivial). It is very important that in 
our case the transfer maps generate the action of an infinite abelian group Z n_1 
so the question of whether two such actions are equivalent or not is a non-trivial 
problem of dynamical systems theory. The fact that the corresponding group is 
abelian indicates that the dynamics is integrable but to make this claim precise 
is not so easy. 

Remark. We should mention that one can introduce another dynamical 
analogue of the transfer matrices, which may be even more natural but leads to 
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the multi- valued maps (correspondences). Indeed consider the map 

F = R 0n Ron-l ■■■Roi ■■ (x ,Xi, ...X n )-> (?/o,l/l, ■■•,Vn) 

and assume that xq — 2/0 (to mimic the trace operation). Generically one 
can express xq as a function of x\, X2, ■ ■ ■ , x n from this relation (in not neces- 
sarily unique way though). Substituting this xq into the rest of the relations 
2/i = Vi(xq, Xi,..., £„), ...,y n = y n (x ,xx, . . . , x n ) we come to the transfer cor- 
respondence v? 1 ' : X n — * X n , which can be considered as analogue of transfer 
matrices Q. The commutativity property (JSj) follows as in the usual case. More 
detail can be found in [52] where the corresponding multi-valued dynamics re- 
lated to Adler maps and discrete KdV equation is investigated. 

Now we are going to present the main constructions and the examples of the 
Yang-Baxter maps coming from the theory of integrable systems. 

Matrix factorisations and Lax representations 

The relation between matrix factorisations and Yang-Baxter property goes back 
at least to A.B. Zamolodchikov's works in the late 70-th on factorised S-matrices 
and related algebras. I would like to mention here also the papers [2,54-58]. 
The first concrete example of the Yang-Baxter map I know which came from 
this construction is the Adler map [34] (see below) . 

Let A(x, Q be a family of matrices depending on the point x £ X and a 
"spectral" parameter ( e C. One should think of X being an algebraic variety 
and A depending polynomially/rationally on £ although the procedure works 
always when the corresponding factorisation problems are uniquely solvable. 
For matrix polynomials usually this is the case once the factorisation of the 
corresponding determinant is fixed (see [53]). 

The following procedure is a version of standard QR-algorithm in linear 
algebra and known to be very useful in the theory of the discrete integrable 
systems, see e.g. [54]. 

Consider the product L — A(y, ()A(x, £), then change the order of the factors 
L — ► L = A(x, C)A(y, £) and refactorise it again: L — A(y, QA{x, Q. Now define 
the map R by the formula 

R(x,y) = (x,y). (14) 

In other words the relation 

A(x, A; QA(y, W = A(y, w QA(x, A; 0, (15) 

is satisfied iff (x, y) = R(X, n)[x, y). 

In that case the matrix A(x 1 £) is called the Lax matrix of the map R, which 
is automatically reversible Yang-Baxter map. Indeed the Yang-Baxter equation 
easily follows from associativity and the reversibility is obvious. 

Once the Lax matrix is known one can introduce the monodromy matrix 

M = A(x n , QA(x n -i, C) • ■ • A{x u C) . 
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Theorem [24]. The transfer maps T^\ 2 = 1,... , n preserve the spec- 
trum of the corresponding monodromy matrix M(x\, . . . , x n ; Q for all £. 

The proof is simple. For the map t{™^ = R\ n Ri n -i ■ ■ ■ R12 °ne can easily see 
that A{x n ,QA{x n _ 1 ,Q...A{x 1 ,Q = A{x u QA{x n , QA{x n - U Q . . . A(x 2 , C), 
where (x±, . . . ,x n ) = T^ n \x\, . . . ,x n ). Since the spectrum of the product of 
the matrices is invariant under the cyclic permutation of the factors, we have 

SpecM(xi, ...,x n ;Q = SpecM(xi, . . .,x n ;Q- 

To prove the same for one should consider the matrix A(x n +i-\, £) . . . A(xi,Q 
which obviously has the same spectrum as M and use the same arguments. 

As a corollary we have that the coefficients of the characteristic polynomial 
X = det(M— XI) arc the integrals of the transfer maps T^ n \ If the dependence of 
A (and therefore M) on A is polynomial then the dynamics can be linearised on 
the Jacobi varieties of the corresponding spectral curves x(C; A) = (sec [54]). 
If these maps are also symplectic then one can use also the discrete version of 
the Liouville theorem [37] to claim their integrability. 

We should mention here that although the procedure is similar to the one 
proposed in [54] for the discrete Lagrangian systems there is one important 
difference: in the construction of Yang-Baxter maps the factors A should depend 
only on one of the variables, while in general scheme they usually depend both 
on x and y. 

Example. Symmetries of the periodic dressing chain and recutting of poly- 
gons [33,34]. 

In [34] Adlcr discovered a remarkable geometric dynamical system, which he 
called the recutting of polygons. It turned out that (after a slight modification) 
it can be described by the following refactorisation procedure related to the 
periodic dressing chain investigated by Shabat and the author [33] . 

Here I = CxC and the matrix A(x), x = (/, (3) has the form: 

^({ /2+ ;- c ) 

One can check that the factorisation procedure described above leads in this 
case to the following birational map R : /2, P2) — * fi, P2) ■ 

7 _ f Pi ~ 02 7 _ f & ~ fk n7 , 

Jl - 12 ~ , , J2 - Jl ~ , , , (17) 

Jl + .12 Jl + J2 

which is the Adler map. Strictly speaking Adler considered the S- version of this 
map: S = PR. The observation that it satisfies the braid relations and thus 
leads to the action of the affine Weyl group was important for him, because this 
was related to the polynomial growth property [37]. In a different context this 
action was discussed in much more details by Noumi and Yamada [50]. 

The transfer dynamics in this example is symplectic (see also below) and 
integrable in the sense of discrete Liouville theorem [37] (see [34]). This follows 
from the theory of the dressing chain [33] which is closely related to the finite-gap 



9 



theory of the KdV equation [35,36]. In fact the transfer dynamics T^' with odd 
n = 2g + 1 can be considered as a discrete version of the Novikov equations [35] 
describing the (/-gap Schrodinger operators (see for details [33,34]). 

It turned out that quite often (in particular, in the previous example) one 
can derive the Lax matrix directly from the Yang-Baxter map. The following 
observation is due to Suris and the author [26] and was inspired by the results 
[45,46] on 3D consistent discrete equations. 

Suppose that on the set X we have an effective action of the linear group 
G = GLn, and that the Yang-Baxter map R(X,/j,) has the following special 
form: 

x = B(y,fi,X)[x], y = A(x, A, fi)[y], (18) 

where A, B : YxCxC — > GLn are some matrix valued functions on X 
depending on parameters A and /i and A[x] denotes the action of A 6 G on 
x <E X. Then the claim is that both A(x, A, () and B T (x, A, £) are Lax matrices 
for R. The claim about B is equivalent to the fact that B(x 1 A, () is a Lax 
matrix for R%x. 

Indeed, look at the values of z 12 produced by the both parts of the Yang- 
Baxter relation |5|): the left-hand side gives z 12 — A{y\, /i, v)A{x 2l A, v)[z\, 
while the right-hand side gives z\ 2 = A(x, A, v)A(y, /i, v)[z\. Now since we as- 
sume that the action of G is effective, we immediately arrive at the relation 
A(x, A, v)A(y, fi, v) = A(y±, u)A(x 2 , A, v), which holds whenever (x 2 ,yi) — 
i?(A, fi)(x,y). This coincides with ([M]) , an arbitrary parameter v playing the 
role of the spectral parameter (. 

Similarly, one could look at the values of x 2 3 produced by the both parts 
of (J5]): the left-hand side gives x 23 = B(z, X)B(y, fi, X)[x], while the right- 
hand side gives x 23 — B(y 3} fj,, \)B(z 2l v, X)[x]. Effectiveness of the action of G 
again implies: B(z,is, X)B(y, fj,, A) = B(j/3, /x, X)B(z2, v, A), whenever (2/3,22) = 
R(ji,v)(y, z). This turns into (fl4|) for the transposed matrices B T (or for the 
inverse matrices -B -1 ); the role of spectral parameter is here played by an arbi- 
trary parameter A. 

In fact the effectiveness of the action can be replaced by the following weaker 
property. Let us call the action of G = GL n on X projective if the action of A on 
X is trivial only if A is a scalar. Similarly we say that A(x, A, Q gives a projective 
Lax representation for the Yang-Baxter map R if the relation (|15| holds up to 
multiplication by a scalar matrix cl, where c may depend on all the variables 
in the relation. One can easily modify the previous arguments to produce the 
integrals for the transfer maps using the projective Lax matrix: all the ratios 
of the eigenvalues of the monodromy matrix are obviously preserved by these 
maps. 

For the projective action the matrices A(x,X,() and B T (x,X, £) give pro- 
jective Lax representations for the corresponding Yang-Baxter maps fj 14|) . In 
practice however for a natural choice of matrices A, B we have actually proper 
Lax representations, as we can see from the example of Adler map (TT7|) . 
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Indeed one can rewrite (|17p as follows 



x 2 + xy - (fx -A) 4/ \ x r i 
■ = A(a;, A, /j) [y] , 



where as before 



A(x,X,C) = ( I 



x 2 + X-Q 



x 



and the group G = GL 2 is acting on CP 1 by Mobius transformations. 

Another interesting example when the Lax representation can be read di- 
rectly from the form of the Yang-Baxter map is related to the theory of solitons, 
in particular matrix solitons. 



Interaction of solitons as a Yang-Baxter map. 

Consider any integrable by the inverse scattering method PDE in 1+1 dimen- 
sions which has multisoliton solutions. Suppose that each soliton has a non- 
trivial internal degrees of freedom described by the set X (which is usually a 
manifold). The soliton interaction has the following remarkable property: the 
result of interaction of three solitons is independent of the order of collisions 
and thus is completely defined by the pairwise interaction. Although the proof 
of this claim in full generality seems to be absent in the literature it is confirmed 
in all known cases (see recent papers [40-42] for the discussion of this). This 
property means that the corresponding soliton interaction map R satisfies the 
Yang-Baxter relation. 

A good example is the matrix KdV equation [59] : 

U t + 3UU X + 3U X U + U xxx =0, 

U is d x d matrix. 

It is easy to check that it has the soliton solution of the form 

U = 2X 2 Psech 2 (Xx - 4A 3 t), 

where P must be a projector: P 2 = P. If P has rank 1 then it should has the 
form P = 1 where £ is a vector in a vector space V of dimension d, 77 is a 

vector from the dual space V* (covector) and bracket (£, 77) means the canonical 
pairing between V and V*. 

The change of the matrix amplitudes P ("polarisations") of two solitons 
with the velocities Ai and A2 after their interaction is described by the following 
map [25,40]: 

R{Xi,X 2 ) : (£d»7i;6)T72) -> (£1,771 ; £2, 772) 
r , 2 A 2 (£1,772) . 2X 2 (&,Vi) , 1Q v 

6=6 + TT -T77 ^2, m=Vl + 7T ~ T77 7^2, (19) 

(Ai - A 2 )(£ 2 ,?72) (Ai - A 2 )(£ 2 ,?72) 
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6=6 + 7T rTTe r£i, ??2 = ?72 + tt ri7i- (20) 

(A 2 -Ai)(Ci,??i) (A 2 - Ai)(Ci,?7i) 

In this example X is the set of projectors P of rank 1 which is the variety 
CjpiV-i x CP^ -1 , and the group G — GL N is acting on the projectors by 
conjugation (which corresponds to the natural action of GL(V) on 7® V"*). 

It is easy to see that the formulas (|19[). (j20|) are of the form ([18]) with the 
matrices 

^o-'+c^-'+c-*' (21) 

(note that again P 2 i = R)- Our results show that the matrix A(P, A, Q gives a 
projective Lax representation for the interaction map. 

One can show that this is actually a genuine Lax representation. We will 
demonstrate this for the class of more general Yang-Baxter maps on Grassman- 
nians following to [25]. 

Let V be an n-dimensional real (or complex) vector space, P : V — » V be 
a projector of rank k: P 2 = P. Any such projector is uniquely determined by 
its kernel K = KerP and image L = ImP, which are two subspaces of V of 
dimensions k and n — k complementary to each other: K © L = V. The space of 
all projectors X of rank k is an open set in the product of two Grassmannians 
G(k, n) x G(n — k, n). 

Consider the following matrix 

2A 

A(P,A;C)=I+— P (22) 
and the related re-factorization relation 

(/ + ^-Pl)(I + 7^2) = (I + + T— T-Pl) (23) 

which we can rewrite in the polynomial form as 

((C-A 1 )I+2A 1 P 1 )((C-A 2 )J+2A 2 P 2 ) = ((C-A 2 )/+2A 2 P 2 )((C-A 1 )I+2AiPi). 

(24) 

The claim is that if Ai ^ ±A 2 it has a unique solution. This follows from 
the general theory of matrix polynomials (see e.g. [53]) but in this case we can 
see this directly. 

Indeed let us compare the kernels of both sides of the relation (|2"4")) when the 
spectral parameter ( — Ai. In the right hand side we obviously have K\ while 
the left hand side gives 

2A 

((A x - A 2 )7 + 2X 2 P 2 )- 1 K 1 = (I + ^-P,)- 1 ^. 

Ai — A 2 



Now we use the following property of the matrix 

A{P, -A; C) = A(P, A; C) _1 (25) 
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to have 

9 A 

k ^^-x^\ 2 P2)Kl - (26) 

Similarly taking the image of both sides of (|24|) at ( = X 2 we will have 

9 A 

L 2 = (/+t V P i) L 2- (27) 

To find if 2 and Li one should take first the inverse of both sides of (f2"3")) . 
use the property (|25[) and then repeat the procedure. This will lead us to the 
formulas: 

2 A 

*a = (J- r-Y^i)^ (28) 

Al T A2 

and 

2 A 

ii = (I+T V p 2)^i- (29) 

Ai — A2 

The formulas (|26I27I28I29[) determine a parameter-dependent Yang-Baxter 
map on the set of projectors. One can easily check that for k = 1 one has the 
formulas (TT91 |2"0")) for two matrix soliton interaction. 

If we supply now our vector space V with the Euclidean (Hermitian) struc- 
ture and consider the self-adjoint projectors P of rank k then the corresponding 
space X will coincide with the Grassmannian G(k, n) : such a projector is com- 
pletely determined by its image L (which is a fc-dimensional subspace in V 
and thus a point in G(k,n)) since the kernel K in this case is the orthogonal 
complement to L. 

The corresponding Yang-Baxter map R on the Grassmannian is determined 
by the formulas 

L 1 = (I+-^-P 2 )L i , (30) 



Ai - A 



2 



9 A 

£ 2 = (/+t V p i)^- (31) 

A2 — Ai 

Yang-Baxter maps related to geometric crystals 

In this section we briefly discuss the examples of Yang-Baxter maps related to 
geometric crystals and tropical combinatorics [13, 14, 19, 22]. We show how one 
can derive the Lax representation for them following to [26] . 

Let X — C n , and define the map R: XxX^XxXby the formulas 



Pj ~ _ Pj- i 



yj = yj- J 5-> j = i,---,", (32) 



where 



fa-\ 



0=1 \fc=l fe=o+l / 
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(in this formula subscripts j + k are taken (mod n)). The map (j3"2"| obviously 
preserves the following subsets invariant: X\ x X^ C X x X, where X\ = 
{(xi, . . . ,x n ) £ X : YYk=i x k — A}- Let us show that the restriction of R to 
X\ x X M can be written in the form (fT5|) . For this, let us embed this set into 



CP" 1 x CP' 



J{x,y) = (z(x),w(y)), z(x) = (1 : Z\ : . . . : 2„_i), w(y) = (w 1 



1), 



Zj = n x fc 

fe=i 



n 

k=j+l 



Then it is easy to see that in coordinates (z,w) the map R is written as 

z = B(y, fi, A) [z] , w = A(x, X, /i) [w] , 

with certain matrices i3, A from G = GL ni where the standard projective action 
of GL n on CP™ -1 is used. A simple calculation shows that the inverse matrices 
are the following matrices well-known in this context [14,19,22]: 



A ^a.A./i) 





Vi 


-1 


. 








\ 







y-2 


-1 . 


















V3 ■ 


















. 


■ ■ Vn-1 


-1 




V 


-A 





. 





Vn 


) 




Xl 





. . 











-1 


#2 


. . 















-1 


x 3 . . 


















. . 


X n — i 







V 








. . 


. -1 


3> n 


) 



(34) 



(35) 



Thus we derived again the Lax matrices directly from the map itself. It is worthy 
to mention that this is closely related to the notion of the structure group Gr of 
the Yang-Baxter map R [14]. It was shown by Etingof in [14] that for the map 
(|32"|) the so-called reduced structure groups G# and can be realised as the 
subgroups of the loop group PGL n (C(X)) generated by the matrix functions 
A^ 1 {x 1 •, A) with x £ X, resp. by B~ 1 (x, ; A) with x £ X . 



Poisson Lie groups and Hamiltonian structure for 
the transfer dynamics 

In this section we describe the Hamiltonian structure for the transfer maps 
following the recent work by Reshetikhin and the author [27] based on some 
general constructions from Weinstein and Xu [3] . The main observation is that 
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if the Lax matrices form symplectic leaves in certain Poisson Lie group then the 
transfer dynamics is Poisson. 

For a good account of the theory of the Poisson Lie group we recommend 
Chapter 5 in the Reiman and Semenov-Tian-Shansky book [60] or original Drin- 
feld paper [61]. 

Recall that Poisson Lie group is a Lie group G with the Poisson structure 
on it such that the multiplication map fj, : G x G — > G is Poisson. Recall also 
that the map between two Poisson manifolds is called Poisson if it preserves the 
Poisson bracket. 

We will need also the notion of Poisson correspondence (or Poisson relation) 
(see e.g. Weinstein [62] and references therein). The correspondence (relation) 
$ from a set X into itself is a multivalued, partially defined map determined by 
its graph T$ , which is a subset of X x X. Now let X — M be a smooth manifold 
with Poisson structure J and the graph T$ is a submanifold of M x M. Let 
us supply the manifold M x M with the Poisson structure J © (—■/)■ Recall 
now that a submanifold N of a Poisson manifold P is called coisotropic if the 
set of functions on P which vanish on N is closed under Poisson bracket. The 
correspondence $ is called Poisson if its graph T$ is coisotropic submanifold of 
M x M. For the usual maps this equivalent to the standard notion of Poisson 
map. If some smoothness condition is fulfilled then the composition of two 
Poisson correspondences is again Poisson. 

For any Poisson Lie group G one can construct the following correspondence 
$g : G x G — > G x G: its graph T$ is the set of (31,52; hi, h 2 ) eGxGxGxG 
such that 

9192 = h 2 hi. (36) 

From the definition of the Poisson Lie group it follows that this correspon- 
dence is Poisson. 

The first simple observation is that this correspondence satisfies the Yang- 
Baxter relation. Note that for the correspondence R: XxX^XxX the 
map R12 is defined by its graph which is a subset of X^ of the form Tr x A, 
where A = (x,x),x S X is the diagonal in X x X. Similarly one can define 
Rij and thus Yang-Baxter correspondence as such R which satisfies the relation 
|l|. Reversibility is a bit more tricky. The problem is with the notion of the 
inverse correspondence the standard formula = Id seems to be too 

restrictive, so we should avoid it. There are two involutions on X x X x X x X: 

7T : (xi, £2; 2/1, 2/2) = (£2, 2:1; 2/2,2/1) 

and 

t : (xi, £2; 2/1, 2/2) = (2/1, 2/2; xi, x 2 ). 

We will call the correspondence R reversible if ttTr — tTr, or equivalently its 
graph is invariant under the involution 

a = ttt : (£1,2:2; 2/1, 2/2) = (2/2, 2/i; X2, xi). 
For the maps this is equivalent to the definition given above. 
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We claim that for any Poisson Lie group G the Poisson correspondence <£>g 
is a reversible Yang-Baxter correspondence. Indeed the Yang-Baxter property 
follows from the associativity of multiplication in G: both left and right sides 
are the correspondences with the same graph given by the relation 

9i929z = h 3 h 2 h 1 . 

Reversibility is obvious. 

In such generality this claim is probably not interesting enough, since we are 
primarily interested in genuine maps but not correspondences. A natural idea 
therefore is to restrict the previous construction to the symplectic leaves in G 
in order to make the re-factorisation relation (|3"6"]) uniquely solvable. 

Suppose that there exist a one-parameter family of cmbeddings of the set 
X as a symplectic leaf in the Poisson Lie group G: (j>\ : X — > G and define the 
correspondence i?(A, /j.) : X x X ^ X x X by the relation 

4>x(x)Mv) = MAM*)- (37) 

Note that the embedding <f>\ induce on X the symplectic structure 

w A = #(w G ), (38) 

where loq is the symplectic structure on the corresponding leaf in G. Let us 
introduce on X x X a symplectic structure as the direct sum oj\ © w^. In a 
similar way we define the symplectic structure fiW on X^ as the direct sum 
of the forms cj, = : 

fi (Ar) =wifflW2ffi...ffiw,v. (39) 

Theorem [3,27]. The correspondence R(\, (ji) defined by ^37]) is a reversible 
Yang-Baxter Poisson correspondence. The related transfer maps {Pp are Poisson 
correspondences with respect to the symplectic structure fl^ N \ 

This construction is already good enough to produce many interesting ex- 
amples of symplectic Yang-Baxter maps, including the interaction of matrix 
solitons [24,25]. To show this we consider the standard Poisson Lie loop group 
G = LGL(n, R) and the symplectic leaves of a special form. 

The standard Poisson structure on the loop group G = LGL(n,H) is given 
by the Sklyanin formula 

{ii(0, ®L 2 (v)} = [-^-,ii(C) ® L 2 (r,)}, (40) 

C - v 

where P\ 2 is the permutation matrix: P(x ® y) — y ® x. 

We will consider the symplectic leaves of the form !/(£) = A + B(. The 
relation (|40| is equivalent to the following set of relations 

{A 1 ,®A 2 } = P 12 (A 1 B 2 -B 1 A 2 ), {Ai,®B 3 } = 0, {S x ,®B 2 } = 0. (41) 

From these relations it follows that the matrix elements of B are in the 
centre of the corresponding Poisson algebra, so B can be fixed. 
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If B is non-degenerate one can assume that B = I is identity matrix. The 
rest of the relations (|4Tj) determine the standard Lie-Poisson structure on the 
matrix group GL(n, R) considered as an open set in the n x n-matrix algebra 
M n which can be identified as a space with gl*{n, R). 

Let Ok be the special symplectic leaf in (^*(n,R) consisting of involutions 
S = 2P — I, where P is a projector of rank k: P 2 = P, S 2 = I. Such a projector 
P is determined by its kernel K = KerP and its image L = ImP, which are 
two subspaces of dimension k and n — k respectively such that K © L = R n . 
Thus X is an open subset of the product of two Grassmannians G(n, k) and 
G(n,n — k). Note that the Grassmannian G(n,n — k) is naturally isomorphic 
to the Grassmannian G(n, k) in the dual space. 

Let X — Ok and define the embedding <fi\ of X into G by the formula 

<j)x(S) = (I + \S = {( - \)I + 2\P. (42) 

It is easy to see that the formulas (|37|) . (|42|) lead to the Yang-Baxter map 
(|26l27l28l29p considered in the previous section in relation with matrix KdV 
equation. In the rank 1 case we have the formulas (|19[) . ([2H)l . 

The corresponding transfer maps (|9]) are symplectic and preserve the 

spectrum of the monodromy matrix 

M(C) = C(Pi, \i,0£(P 2 , A 2 , • ■ ■ £(Pn,\n,(), (43) 

where C(P, A, C) = (C — A)/ -I- 2XP. The coefficients of the characteristic polyno- 
mial x(Ci v) = dct(£ — r/I) give the set of integrals for transfer dynamics, which 
are in involution with respect to the symplectic structure Q^ N > on X^ N \ 

We believe that this implies integrability of the transfer dynamics in Liouville 
sense, but the proof of completeness is still to be found. As usual the dynamics 
is linear on the Jacobi variety of the spectral curve x(Ci 7 l) = 0- 

Let us consider now the symplectic leaves in LGL(2,R) of the form L = 
A + B( with degenerate 

1 




B = 
Writing A as 



A = 



a b 
c d 



we have the following Poisson brackets for the matrix elements from the relation 
CP: 

{a,b} — a, {a,d} — c, {a, c} = 0, {6,c} = 0, {b,d} — d, {c, d} = 0. 

It is easy to check that the centre of this Poisson algebra is generated by two 
functions: C\ = c and C2 = ad — be. If we choose C\ = c = 1 and C2 —ad — bc — 
A we come to the symplectic leaf 0(A) consisting of the matrices of the form 

MCA)=( 1 ), (44) 
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where {a, d} = 1. 

This is not quite Lax matrix (fP6|) for the Adler map yet. To have it one 
should apply the Hamiltonian reduction to the product of the corresponding 
orbits 0(A) (see [27] for the details). This leads to the Poisson structure for 
the transfer dynamics for Adler maps, which is known from the theory of the 
periodic dressing chain [33,34] (with an odd period N). 



Classification of Yang-Baxter maps: quadrirational 

case 

The classification of the Yang-Baxter maps is probably the most interesting 
but hopelessly difficult problem even in the simplest case of X = CP 1 . A 
natural class where the complete classification seems to be realistic is given by 
the quadrirational maps introduced and investigated by Adler, Bobenko and 
Suris [43]. 

The map R : (x, y) — > (u, v) is called quadrirational if both R and so called 
companion map R : (x, v) — > (y, u) are birational maps of CP 1 x CP 1 into itself. 
All such maps have the form 

a(y)x + b(y) = A(x)y + B(x) 
• U c(y)x + d(y)> V C(x)y + D(x)' 1 ) 

where a(y), . . . , d(y) and A(x), . . . , D(x) are polynomials of degree at most 2. 
There are three subclasses of such maps, denoted as [1:1], [1:2] and [2:2], cor- 
responding to the highest degrees of the coefficients of both fractions in (|45|) . 
In [43] the classification of such maps was done modulo group (Mob) 4 , where 
Mob is the group of projective transformations of CP 1 : 



In the richest class [2:2] it looks as follows. 

Theorem [43]. Any quadrirational map on CP 1 x CP 1 of subclass [2:2] is 
equivalent, under some change of variables from (Ai '6b) i , to exactly one of the 



following five 


maps 






u = ayP, 


v = 


(3xP, 


p 


u — — P, 

a 


V = 


-P 

P ' 


p 


y v 

u = — P, 

a 


v = 


-P 


p 


u = yP, 


V = 


xP, 


p 


u = y + P, 


v = 


x + P, 


p 



(1 - (3)x + (3 - a + (a - l)y 
(3(1 - a)x + (a - (3)yx + a((3 - l)y ' 
ax — (3y + (3 — a 

x-y 
ax — (3y 



x-y 
(3-a 



x-y 

A 

V 



(Fi) 
(Fu) 
(Fm) 
(Fw) 
(Fv) 
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with some suitable constants a, (3. 

It is remarkable that all these maps can be described by the following natural 
geometric construction [43]. 

Consider a pair of nondegenerate conies Q\, Qi on the plane CP 2 , which 
are rational curves isomorphic to CP 1 . Take X £ Q\, Y G Q2, and let I = XY 
be the line through X,Y. Generically the line I intersects Q\ and Q2 at two 
points U and V (see Fig. [BJ. This defines a map K : {X, Y) i-> (U, V). 




Figure 6: A quadrirational map on a pair of conies 

There are five possible types I — V of intersection of two conies (see [63]): 
I: four simple intersection points; 

II: two simple intersection points and one point of tangency; 
III: two points of tangency; 

IV: one simple intersection point and one point of the second order tangency; 
V: one point of the third order tangency. 

In suitable rational parametrisations of the conies we have exactly the five maps 
listed in the theorem above. 

It is remarkable that all five representative maps from the last theorem 
satisfy the Yang-Baxter relation (see [43] , where a nice geometric interpretation 
of this relation is also given). This surprising fact, which must be related to 
the classification approach of [43] based on the analysis of the singularities, 
does not imply however that all quadrirational maps are Yang-Baxter (as one 
might conclude from [43]) since the action of (Aidb) 4 in general destroy the 
Yang-Baxter property. 

For example, if we change in the fifth map Fy the variables x — > —x, y — » — y 
(preserving u, v) we come to the quadrirational map u = —y—P, v = — x~ P, P = 
^rr--, which does not satisfy the Yang-Baxter relation. However, if we change 
u — » — ii, y — » —y we come to the Adler map (fTTj) . which is a Yang-Baxter 



19 



map. Note that neither of these changes is a proper conjugation of maps of 
CP 1 x CP 1 into itself (in particular, Fy is not conjugated to Adler map). 

A question remains therefore to classify all quadrirational Yang-Baxter maps 
modulo diagonal action of Mob, which seems to be a natural equivalence within 
this class. The previous example shows that in the same (Mob) -orbit there are 
several Yang-Baxter representatives, non-equivalent under A4o6-conjugation. 
The answer to this question can be probably extracted from the results of the 
paper [43]. 

Let us mention in this relation that a complete classification of intcgrablc 
(3D-consistent) discrete equations on quad-graphs is known due to the same 
authors [44]. Since these two problems are closely related (see [44,65]) a very 
interesting question is what could be a Yang-Baxter analogue of the "master 
equation" Q4 in their classification (see [64] for the discussion of this very impor- 
tant case). Probably one should look for a wider class of Yang-Baxter maps (like 
correspondences or Cremona maps). A class of matrix factorisations suggested 
by Odesskii [58] could give one of the possible answers to this question. 
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